The Intrinsic Difficulties of Constructing Strongly Correlated States of Lattice 
Quantum Gases by Connecting Up Pre-engineered Isolated Atomic Clusters 
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Suppose one engineers an artificial antiferromagnet on an array of isolated wells, and then increases 
the tunneling between wells slowly, will the system finally become an equilibrium antiferromagnet? 
Here, we show that due to the intrinsic non-adiabaticity at the start of this process, and that the 
atoms in the initial state in different wells are completely uncorrelated, the final equilibrium state 
will have a temperature Tf far above the Neel temperature. Constructing other strongly correlated 
states (with characteristic energy per particle E* comparable to the hopping or the virtual hopping 
scale) with the same method will suffer the same problem, i.e. Tf >> E* . 



At present, there are intense experimental efforts to 
manufacture strongly correlated states with cold atoms in 
optical lattices. The success of these attempts will allow 
one to study many challenging models in condensed mat- 
ter physics. These efforts, however, are exceedingly chal- 
lenging. The reason is that the energy scales associated 
with these correlated states are very small, corresponding 
to temperatures as low as 10 _12 K and beyond[l|. Reach- 
ing these unprecedentedly low temperatures is a great 
achievement in itself. 

With the problem of cooling looming at the horizon, 
other ideas of achieving strongly correlated states in lat- 
tice quantum gases begin to emerge. One popular idea 
is to engineer the desired many-body state directly, do- 
ing away with cooling completely. Within this scheme is 
also the idea of connecting up isolated clusters of atoms, 
previously engineered into specific quantum states. The 
hope is that by creating the right kind of clusters, one 
might achieve the desired equilibrium many-body states 
by connecting them up. For brevity, we shall referred 
this scheme as "direct construction method" . 

For example, to manufacture a resonant valance bond 
(RVB) liquid or RVB solid, one might start with an ar- 
ray of isolated double wells, each of which contains a 
singlet pair of spin- 1/2 fermions, as in the experiment 
of S. Trotzky et.al.[2j One then adjusts the barriers be- 
tween neighboring double wells so that at the end one 
has a cubic lattice with one fermion per site. The hope is 
that if this process is slow enough, the initial state may 
evolve adiabatically into a RVB liquid or a RVB solid. 
Another example is to manufacture an antiferromagnet 
using spin- 1/2 fermions in an optical lattice. One could 
first engineer an "artificial" antiferromagnet consisting 
of alternating up and down spins in a cubic array of iso- 
lated wells Q. It is tempting to think that by reducing 
the barriers between wells sufficiently slowly, the system 
will turn into an equilibrium antiferromagnet in a cubic 
lattice. After all, the initial state is already an "antifer- 
romagnet" with essentially zero entropy. 

The purpose of this paper is to point out an intrinsic 
difficulty of this direct construction scheme. We shall see 



that the process of connecting up the isolated clusters 
is intrinsically non-adiabatic. This non-adiabaticity, and 
the lack of correlations between atoms in different wells in 
the initial state, will lead to an energy difference ~ NE* 
between the initial and the final equilibrium state, where 
N is the number of particles, and E* is the characteris- 
tic energy per particles of the strongly correlated state. 
This excess energy will lead particle fluctuations at the 
surface, so much so that the temperature Tf of the equi- 
librium state is far above the characteristic energy scale 
E* , melting away the correlation one set out to achieve. 
In the case of anti-ferromagnet, E* is the Neel tempera- 
ture Tn, and Tf >> Tjy. While our proof does not ex- 
clude the possible success of other engineering schemes, it 
illustrates the kind of problems that must be faced, and 
the stringent conditions that must be satisfied to reach 
an equilibrium strongly correlated state. As of now, solv- 
ing the problem of cooling, (or more correctly, reducing 
the entropy of the system), remains essential in achiev- 
ing strongly correlated states of cold atoms in optical 
lattices Q . 

I. The initial state and the process of equilibra- 
tion: Our initial state is an artificial antiferromagnet 
on a cubic lattice of isolated sites, 



i* >= LR T II a >>< 



(1) 



where La and Lb are the two interpenetrating sub- 
lattices of a cubic lattice. The superscript ' means only 
sites within a radius R from the center are occupied. The 
density profile therefore has a sharp edge. 

The process of connecting up the isolated sites can be 
described by the following (time dependent) hamiltonian 
Hj = HjQJ) + V + T, where Hj(U) is a 3D Hubbard 
model describing the internal energy of the a two com- 
ponent Fermi gas in an optical lattice, 
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and V is the potential energy of a harmonic trap 

V = V r (u))fi r , V r (u) = ^Mcu 2 r 2 (3) 

r 

Here, aJ CT is the creation operator of a fermion with 
pseudo-spin a at site r of a cubic optical lattice, (a =|, j), 
h ra = al a a ra is the number of fermions at site r with spin 
a, n r — h r <f + n T i; U is the on-site repulsion, J is the tun- 
neling integral between nearby sites, and V T — ^Mu 2 r 2 is 
an harmonic potential with frequency lo. The tunneling 
integral J is time dependent. It grows from to a finite 
value over a period of time as the lattice height is de- 
creased from infinity to a finite value. J = corresponds 
to isolated clusters. 

The term Y is the random perturbation due the en- 
vironment. It can be caused by the tiny fluctuations 
in the laser field, or the noise in the current producing 
the magnetic trap, etc. Such perturbations are typically 
small, and are certainly much weaker than, say potential 
energy, i.e. (r 2 ) 1 / 2 << (V). However, they can cause 
de-coherence in quantum evolution. Such de-coherence 
effect are particularly important at the beginning of the 
connection process, as the virtual hopping scale J 2 /U is 
vanishingly small and will be dominated by any small 
perturbation Y. Because of the strong de-coherence at 
the beginning, the system at a later time cannot be de- 
scribed as simply a quantum mechanical evolution of ini- 
tial state l^o) under Hj(U). Rather, one should discuss 
the properties of the system in terms of density matrix, 
and how a density matrix associated with the initial state 
l^o) relaxes to equilibrium. 

Before proceeding, let us quantify the properties of 
jv^o). The number density of Eq.JT]) can be written as 

r4 o) = Oifio - V r ), Mo = Mlo 2 R 2 J2 (4) 

where \l is the "initial" chemical potential, related 
to total particle number as N = J2 r n r°\ or N = 
(27r/3)[/x /(Mw 2 a 2 /2)] 3 / 2 , where a is the lattice spacing. 
To eliminate double occupancy, we also have /i < U . 

II. Formulation of the problem: The evolution of 
the density matrix will be controlled by : 

(A) J depends exponentially on the barrier height B, 
J oc e~ B . As B reduces from infinity to a finite value, 
J remains exceedingly small for a while. During this pe- 
riod, the tunneling time remains so long that the system 
can not be in equilibrium and the process will not be 
adiabatic.H 

(B) Equilibrium only sets in when J becomes sufficiently 
large, which only occurs at a certain time r after the 
infinite barrier was lowered. Because J is changing at a 
finite rate during this period, energy is being deposited 
into the system. The determination of this added energy 
is a complicated problem. 

(C) In practice, the time r to reach equilibrium must be 
shorter than the lifetime of the sample, which is limited 



by three body collisions and other factors. This imposes 
a lower bound on the swept rate of J, and forces one to 
face the problem mentioned in (B). 

While both (A) and (B) are intrinsic, (C) comes 
about because of the specific type of atoms (i.e. the al- 
kalis) used in current experiments. To simplify matters, 
let us assume that we have atoms with infinite lifetime, 
and conduct the following thought experiment: 

(i) We lower the barrier between wells from infinity to a 
very large but finite value J over a time interval t Q . The 
tunneling parameter J is sufficiently small so that the 
tunneling time h/ J is much longer that the switching tim 
t , i.e. t a « h/J. On the other hand, J is sufficiently 
large so that even the smallest energy parameter of the 
system, J 2 /U , will from now on dominates over the noise 
term Y, i.e. J 2 /U >> 7. We therefore have following 
sequence of time scales 

t « h/j « h/(j 2 /u) « n/ 7 . (5) 

(ii) After that, we wait for the system to reach equilib- 
rium. Due to the low tunneling rate, this process will take 
a very long time. However, since the atoms have infinite 
lifetime, equilibrium will set in eventually. What we are 
interested in are the temperature T and the entropy S of 
the final equilibrium state. 

Since process (i) takes place much faster than the tun- 
neling time, the quantum state at the end of this process 
becomes j*,*) = \$> ) + 0(e), where e = -ft /h « 1, (see 
eq.flU). As a result, the internal energy and potential 
energy are E\ nt — (Hj(U))^ i and Vi = (V}* 4 are given 
by the value of \^> ) to the zeroth order in e, i.e. the total 
energy at the end of process (i) is (or order 0(e)), 

Ei = Ef 1 + V t - Vr(u)n r o) = 0. (6) 

r 

Next, we note that energy is conserved in process (ii). 
We then have E t = E f , i.e. EJ Lt + V f = Ej nt + V t , or 

V f -V.^Y. W(n, ~ 4 o) ) - El nt - Ef, (7) 

r 

where n r and Ef are the density and internal energy 
of the final equilibrium state, which are functions of the 
final temperature T. Eq.j7|) shows that the energy differ- 
ence between the initial state and the final equilibrium 
state will lead to particle excitation at the surface, which 
is the origin of entropy production. Eq.([7]) allows us to 
determine T of the final state once the temperature de- 
pendences of E mt and n T are known. Finally, we note 
that for high barriers and no double occupancy, we have 

J 2 /U « J «fi a < U. (8) 

III. The energy of the final equilibrium state : 

We shall assume the final temperature T is less than U, 
for otherwise the system will not be magnetically ordered. 
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For large barrier heights, U >> J, the term J in eq.([2]) 
can be treated as a perturbation. To the zero order in J, 
we have 

2 e (M-V r )/T + 2e [2( i u-V r )-C/]/T 

nr ( T ) = i + 2e(^-^)/ T + e [ 2 (M-^)-c/]/T ' ( 9 ) 

where \i is chemical potential of the final equilibrium 
state, determined by N = ^ r n r . Since the chemical po- 
tential of initial state satisfies fi < U, we expect n < U, 
i.e. double occupancy remained disfavored. In this case, 
(which we shall verify later), we have 



n r (T) 



l e (%-M)/T + 1 ' 



N 



(10) 



Eq. (fT0|) is essentially a step function with size R and a 
step width AR, where 



H = Mcu 2 R 2 /2, AR = T/(Moj 2 R). 



(11) 



The system can then be divided into following regions: 

Region 1 : r < R- AR/2, n T = 1 (12) 
Region2 : R - AR/2 < r < oo, < n r < 1 (13) 

Let Ei and E? be the total internal energy in region 1 
and 2, and AN be the number particles in region-2. The 
internal energy E\ nt of the final state is 



E) nt = Ex 



E 2 , Ex = e(T)(N — AN), 



E 2 = e^AN, 
(14) 

where e(T) is the energy per particle of a bulk equilibrium 
anti-ferromagnet at temperature T, and is the average 
energy per particle inside the mobile layer. 

To find the energies e and €2, we note that in the limit 
of U >> J, the Hubbard hamiltonian Hj(U)( is reduced 
to the tj model H t j =f + Hj®, which is a sum of 
a Heisenberg interaction Ttj = ^3J2{r r') ■ S r ', J = 
J 2 /U, and a correlated hopping term T = —Ja\ iy a r , (71 
where S r is the psuedo-spin 1/2 operator of the fermion 
at site r, and aJ CT is a "correlated" creation operator 
within the space of no double occupancy. 

In region-1, the system has one fermion per site. Hence 
T = 0, and H t j reduces to the 3D anti- ferromagnetic 
Heisenberg hamiltonian Hj. According to the latest 
estimate [7[, the Neel temperature is Tjv = ^3 with 7 = 
0.944. On dimensional grounds, e(T) = -[a(T/T N )]T N , 
where a is dimensionless function of T/Tn. Standard 
spin wave calculation [8] and high temperature series ex- 
pansion give: 

3 58 T — > 

a ^ 7 - 2 (3/4) 2 7VT Wh6n T»T N - (15) 



(Ti-j) < 0. Hence, we have 62 < 0. In addition, since 
{Hj) ~ {J 2 /U)x, and (T) - J(l - x), where x is the 
average number of fermion per site in region-2. Hence 
we have leil/J ~ (1 — x ) + (J/U)x ~ 1. Combining 
eq.© and l]14p. we have 

5^ V r (w)(n r - n( o) ) = aTAr(iV - AiV) + |e5|AiV. (16) 

r 

Since Tn ~ J 2 /?7 << J, and (eil/J ~ 1, we have the 
following condition on the temperature of the final state, 

JN > Vr(oj)(n r - n^) > aT N (N - AN). (17) 



(18) 



Note that n T also satisfies 



IV. The temperature of the final state: Eq. (jT8j) 
and eq. (|16p determine T and (i of the final state in terms 
of the chemical potential /i Q of the initial state, should the 
temperature dependence of a and I2 are known. Since 
eq. ([T0| reduces to a step function 6(fx — V r ) as T — > 0, one 
can apply Sommerfeld expansion to evaluate the sums in 
eq. tTrg)) and eq. p^)) . We obtain (see Appendix) 



/j, = Ho 1 - B- 



T 

/'-o 




3A\ (T 



(19) 



(20) 



where B = 0.69 and A = 1.64. If AN « N, (which we 
verify later), Eq.® and ((TTJ) imply that 



(21) 



From the temperature dependence of a in ea. (|15|l . and 
the relations in eq.©, Eq. ([2T|) can be written as 



1 1 — 
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where the (Tn/^o) 1 / 2 and (Tjv/mo) 1 ^ 3 dependence cor- 
respond to the low and high temperature behavior of a, 
and we have omitted proportional constants of order 1. 
The last inequality of Eq. ([22]) shows that 



T > 



Ts [ 

-In 



1/2 



or T, 



N 



Mo 
T N 



2/3' 



»T N (23) 



In region-2, the mobile layer, the sites arc cither empty 
or singly occupied. In this case, (T) < 0. Furthermore, 
due to short range anti-ferromagnetic order, we also have 



The last inequality follows from Eq.® and the fact 
that T N = 0.944(J 2 /C/)|1, which implies fi a /T N 



Ho/{J /U) » 1. The final equilibrium state is therefore 
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not spin ordered. It is a paramagnetic (rather than anti- 
ferromagnet) with entropy per particle ~ ln2 or larger. 
From Eq.(|22]). the assumptions fi < U and AiV « N, 
from which Eq. (f2"Tj) is derived, are also easily verified Q- 

V. The condition for "direct construction" to 
work: From our analysis, one sees that the severe heat- 
ing shown above is caused by the energy difference be- 
tween the initial state and the final equilibrium state, 
which is of order NT^. This excess energy produces par- 
ticle excitations at the surface, which generates entropy. 
This energy difference is inherent in the direct construc- 
tion scheme. The very fact that one starts with isolated 
clusters means all the correlations between clusters in the 
final equilibrium state are all missing in the initial state, 
which produces an energy difference proportional to the 
number of clusters, which is of order N . 

One might think that if the traps are tighten imme- 
diately after the barriers between neighboring wells are 
lowered, one will reduce particle fluctuation and hence 
the final temperature. Mathematically, it means chang- 
ing Eq.© to '£ ir (V I (u')n r -V I (u)ni o) ) = -Ef\ J > u, 
which will certainly change our conclusions. This, how- 
ever, does not work because tightening the trap in the 
non-adiabatic regime will simply change the initial state 
to a new one that sees a tighter harmonic trap. In others 



words, Eq.([7]) will become J2 r U-( w ')( n r 



» 



rpint 



rather than the expression mentioned above. We are then 
back to the previous situation, (i.e. Eq.(|7|) with u) re- 
placed by to'), from which the severe heating follows. 

While our discussions are for antiferromagnets, our 
considerations also apply to similar schemes for other 
strongly correlated states characterized by an energy per 
particle E* , which is the energy (or temperature) scale 
below which the strong correlations between particles 
emerge. Since the energy of particle fluctuation at the 
surface is ~ (T 2 /fi)N, in order for final temperature T 
to be less than E* , T = (E*,( < 1, one needs 



C Z E* 



CE* 



1 2 = E* 2 (24) 

Since /i >> E* in most cases, this condition is very strin- 
gent indeed. 

Appendix: Derivation of Eq. (fl9|) and ([20]) Defin- 
ing £ = V r = MwV/2, we have J2 r = J-^ 
where D{£) = a - 3 dr/d£ = C\[£ for £ > and D(£) = 
for £ < 0, C = it /(Mu) 2 a 2 /2) 3 / 2 and a is the lattice spac- 
ing. Eq. (fl8|) and IT6|) can be written as 

d£D(£)n(£) = / d£D(£) = (2/3)n D(/j, ) 

- OO J — OO 

(25) 



d£D(£)£n(£)- / d£ D{£)£ (26) 



V f - Vi 



where Vf — Vi = J2r V T {uj){n T — n^), and n(£) — 
(ie( £ ~^/ T + l) -1 . Since n(£) is close to but not ex- 
actly the Fermi function, we need to generalize the usual 



Sommerfeld expansion. For any function H{£), we have 
d£H{£)n(£)= / d£H{£) + H' {£*)T 2 A + 0{T 3 ) 

-OO J — OO 

(27) 

where £* is defined as 

£ * = L £ (- d £) d£= » +BT (28) 

with B = /"„ V M= IZ, Tl^dx = 0.69, 



dx 



1.64. Applying Eq. ([27) to Eq. lpgjl . we have 

{£* -n )D(£*) + D'(£*)T 2 A = 0. (29) 
Together with Eq. ((28|) . we have 



Mo - A 



D'(£*) , 

) / T 2 = a + BT. 
D{£*) M 



(30) 



This gives eg. (|19|) upon iteration. Likewise, we obtain 
from Eq.flHl) and (|27|) 



(31) 



V f -V 



f D{£)£ + [D(£*)£*]'T 2 A. 



Using Eq.©, and (30]), we have V f - V = 

D{£*)T 2 A = (3/2)AN(T 2 /fj, ) + 0{T 3 /fi 2 o ), which is 
Eq.pl). 
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